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Abstract 

It is suggested that the dynamics of liquid water has a component 
consisting of 0~ 2z (oxygen) anions and H +z (hydrogen) cations, where 
z is a (small) reduced effective electron charge. Such a model may 
apply to other similar liquids. The eigenmodes of density oscillations 
are derived for such a two-species ionic plasma, included the sound 
waves, and the dielectric function is calculated. The plasmons may 
contribute to the elementary excitations in a model introduced recently 
for the thermodynamics of liquids. It is shown that the sound anomaly 
in water can be understood on the basis of this model. The results 
are generalized to an asymmetric short-range interaction between the 
ionic species as well as to a multi-component plasma, and the structure 
factor is calculated. 

Introduction. As simple as it may appear, water is still a complex liquid 
involving various interactions as well as kinematic and dynamic correlations. 
It is widely agreed that the water molecule in liquid water preserves to some 
extent its integrity, especially the directionality of the sp 3 -oxygen orbitals, 
though it may be affected substantially by hydrogen bondsQ As such, it is 
conceived that water has a molecular electric moment, an intrinsic polariz- 
ability and hindered rotations (librations) which may affect its orientational 



1 L. Pauling, General Chemistry, Dover, NY (1982); Water: A Comprehensive Treatise, 
ed. by F. Franks, Plenum, NY (1972). 
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polarizability. We examine herein another possible component of the dynam- 
ics of the liquid water, as resulting from the dissociation of water molecule. 

Water molecule H 2 has two H — O (hydrogen-oxygen) bonds which make 
an angle of cca 109° in accordance with the tetragonal symmetry of the four 
hybridized sp 3 -oxygen orbitals. The "spherical" diameter of water molecule 
is approximately 2.75A and the inter-molecular spacing in liquid water un- 
der normal conditions is a ~ 3A. This suggests that water molecule in liquid 
water, while preserving the directionality of the oxygen electronic orbitals, 
might be dissociated to a great extent. Dissociation models which assume 
OH~ — H + or OH~ — H 3 + pairs are known for water. This indicates 
a certain mobility of hydrogens (and oxygens). We analyze herein the hy- 
pothesis that water may consist of 0~ 2z anions of mass M = IQamu and 
density n and H +z cations (protons) of mass m = lamu and density 2n, 
where z is a small reduced effective electron charge (the atomic mass unit is 
lamu ~ 1.7 x lCT 24 g.). We shall see that such a hypothesis adds another di- 
mension to the dynamics of water. Such a model may apply to other similar 
liquids. 

Due to their large mass the ions have a classical dynamics. Herein, we limit 
ourselves to considering the ions motion in water under the action of the 
Coulomb potentials (poo — 4:Z 2 e 2 /r, (Phh = z 2 e 2 /r and ipon = —2z 2 e 2 /r, 
where — e (~ —4.8 x 10~ 10 esu) is the electron charge and r denotes the 
distance between the ions. For stability, it is necessary also to introduce a 
short-range repulsive (hard-core) potential x@ It is shown that in the limit 
z — > water may exhibit an anomalous sound-like mode beside both the 
ordinary (hydrodynamic) one and the non-equilibrium sound-like excitations 
governed by short-range interactions. We compute the density oscillations 
for this model, the dielectric function, the structure factor, and extend the 
model to a multicomponent plasma, including an asymmetric short-range 
interaction between ion species. 

Plasmons in a jellium model. Let us consider one species of charged par- 
ticles, with charge — ze, continuously distributed with density n in a neutral- 
ising rigid continuous background of positive charge. This is the well-known 
jellium modelj^l The Coulomb interaction reads 

U = ^J drdr'cp(r - r')Sn(r)5n(r') , (1) 



2 See also in this respect E. Teller, Revs. Mod. Phys. 34 627 (1962); E. H. Lieb and B. 
Simon, Phys. Rev. Lett. 31 681 (1973); Adv. Math. 23 22 (1977); L. Spruch, Revs. Mod. 
Phys. 63 151 (1991). As it is well-known, a classical plasma with Coulomb interaction 
only is unstable. 

3 See, for instance, D. Pines, Elementary Excitations in Solids, Benjamin, NY (1963). 
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where 5n(r) denotes a small disturbance of density (which preserves the 
global neutrality). We introduce the Fourier representation 

Sn{r) = -^^5n(q)e* qr , 5n{q) = I dr<5n(r) e - iqr , (2) 
v N q v N J 

where N = nV is the total number of particles in volume V. Similarly, 

= ^E^) eiqr > ^ = I dMr)e~^ , (3) 
q J 

where y?(q) = 4irz 2 e 2 /q 2 is the Fourier transform of the Coulomb potential 
(interaction). The Coulomb interaction given by (JTJ) becomes 

tf=^I>(<z)Mq)<K-q) ( 4 ) 

in z — ' 
q 

(where the q = 0-term is excluded by the positive background). 

The small variations 8n(r) in density can be represented as 5n = —ndivu, 
where u is a displacement vectorQ We emphasize that such a representation 
holds for qu(r) <C 1. It follows Sn(q) = — mqu(q), and one can see that the 
Coulomb interaction involves only longitudinal components of the displace- 
ment vector u(q) along the wavevector q. Therefore, we may write u(q) = 
(q/g)u(q), with 5n*(— q) = 8n(q), u*(— q)=u(q) and u*(— q) = — u(q). The 
Coulomb interaction (j3J) becomes 

u = ~ ?V(9)w(q)w(-q) • (5) 

q 

The kinetic energy associated with the coordinates u(q) is given by 

T = — J drnmu 2 = — m u(q)u(— q) , (6) 
2 J 2 q 

where m denotes the particle mass. The equations of motion obtained from 
the Lagrange function L = T — U are 

mu(q) + nq 2 (f(q)u(q) = , (7) 

which leads to the well-known plasma oscillations with frequency given by 

= 4irnz 2 e 2 jm. 



'M. Apostol, Electron Liquid, apoma, MB (2000). 
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Plasma oscillations with two species of ions. We apply the above model 
to the two species of ions 0~ 2z and H +z . The change in density is associated 
with a displacement vector v in the former and a displacement vector u 
in the latter. First we note that the Fourier transforms of the Coulomb 
potentials are given by <foo = 4y?(g), tp H H = <f(q) and ip H = -2(p(q), 
where ip(q) = An z 2 e 2 / q 2 . Therefore, the interactions can be written as 

Uoo = -f E q q 2 IMq) + x{q)} v(q)v(-q) , 
u HH = -2nE q 9 2 + x(g)MqM-q) , (8) 
u H = nEq^Pvl?) - xMoM-q) , 

where n = N/V is the density of water molecules and the Fourier transform 
X of a hard-core potential has been introduced (the same for both species). 
The kinetic energy is given by 

T = ~\ M 2 - m Yl ' • ( 9 ) 

q q 
and the equations of motion read 

mil + 2nq 2 ((p + x) u — nq 2 (2ip — x)v = 

Mil + nq 2 (4(f + x)v — 2nq 2 (2ip — x)u = , 
where we have dropped out the argument q. 



(10) 



The solutions of these equations can be obtained straightforwardly. In the 
long wavelength limit q — ► there are two branches of eigenfrequencies, one 
given by 

< = " (11) 

corresponding to plasma oscillations and another given by 

2 $ n X 2 2 2 / 1n \ 

u. = — — - — q = v*q (12) 

s M + 2m s v ; 

corresponding to sound-like waves propagating with velocity v s given by (|T2l) . 
H = 2mM/ (2m + M) is the reduced mass. The plasma oscillations are 
associated with antiphase oscillations of the relative coordinate {2mu + Mv = 
0), while the sound waves are associated with in-phase oscillations of the 
center-of-mass coordinate (u — v — 0). 
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Polarization. An external electric field arising from a potential 0(r) gives 
an additional energy 

Ui = qi / drcj)(r)5ni(r) = -i(n i q i /n) 1 ^q<p(c\)u i (-q) , (13) 
J q 

for a species of ions labelled by z, with electric charge % and density rij. We 
apply this formula to the two-species ionic plasma, and get 

U H = -2ize #(q)«(-q) , U = 2ize qM<l) v (-<l) ( 14 ) 



Adding these two terms to the lagrangian, the equations of motion given by 
(fTUl) become 

mil + 2nq 2 ((p + x)u — nq 2 (2ip — \)v = —izeqtfi 

(15) 

Mv + nq 2 (4(p + x) v — 2nq 2 (2ip — x)u = 2izeq<j) , 

where we have dropped out the argument q. This is a system of coupled 
harmonic oscillators under the action of an external force. In the limit of 
long wavelengths its solutions are given by 



2_ 2m 2 
m ' (u) 2 — ui 2 )(uj 2 — w 2 ) ' 



, ,2 2M , ,2 



(16) 



On the other hand, equation riidivui = —5rii is in fact the Maxwell equation 
divEi = Arrqidrii, where the electric field is given by Ej = —AimqiUi. We have 
therefore the internal electric fields E u = Snnzeu and E v = 8imzev. The 
polarization P = —{E u + E v )/4ir is given by 

2 

P(q) = 2n^e [ U (q) - v(q)] = -^0(q)^^ . (17) 

47r — Up 

The external field is related to the external potential through D(q) = — iq<p(q) 
and the dielectric function e is given by D = eE = s(D + Ew), where 
E int = E u + E v is the internal field. We get the dielectric functional 

e = l- uj 2 /uo 2 , (18) 



3 We disregard here the intrinsic and orientational polarizabilities. 



5 



as expected. As it is well-known, its zero gives the longitudinal mode of 
plasma oscillations. 

The u p in the nominator of equation (fl"8l) defines also the plasma edge: for 
frequencies lower than u p the electromagnetic waves are absorbed (the refrac- 
tive index is given by n 2 = e). It is well-known that water exhibits indeed 
a strong absorption in the gigahertz-terrahertz regionj^ On the other hand, 
neutron scattering on heavy waterjll as well as inelastic X-ray scatteringjj 
revealed the existence of a dispersionless mode ~ 4 — bmeV (~ 10 13 s _1 ) 
in the structure factor, which may be taken tentatively as the u; p -plasmonic 
mode given by equation (fTTl l. Making use of this equation we get u p ~ 
3 x 10 14 zs~ 1 (n = 1/a 3 , a = 3 A), so we may estimate the reduced effective 
charge 2~3x 10~ 2 . 

Dielectric function. The dielectric function given by equation (TTSj) has 
a singularity for u = 0, as arising from the exact cancellation in the static 
limit of the external field by the internal field. It is plausible to assume 
that residual polarization fields are still present in this static limit, like, for 
instance, the intrinsic polarizability. In this case, equation (fLSl) is modified, 
and the dielectric function is of the type 

where u is a plasma frequency associated with the intrinsic, molecular po- 
larizability^] As such, it is a very high frequency, and equation (fT9l) gives 
a small, negative contribution to the dielectric function in the static limit 
(w-0). 



6 See, for instance, K. H. Tsai and T.-M. Wu, Chem. Phys. Lett. 417 390 (2005); A. 
Padro and J. Marti, J. Chem. Phys. 118 452 (2003); K. N. Woods and H. Wiedemann, 
Chem. Phys. Lett. 393 159 (2004). 

7 F. J. Bermejo, M. Alvarez, S. M. Bennington and R. Vallauri, Phys. Rev. E51 2250 
(1995); C. Petrillo, F. Sacchetti, B. Dorner and J.-B. Suck, Phys. Rev. E62 3611 (2000). 

8 F. Sette, G. Ruocco, M. Krisch, C. Masciovecchio, R. Verbeni and U. Bergmann, Phys. 
Rev. Lett. 77 83 (1996). 

9 A static field D produces an electric dipole p = q e x, where q e is the electric charge and 
X is a small displacement subjected to the equation of motion m e x + m e ujpX — q e D, where 
m e is the mass of the electronic cloud. According to the plasma model suggested here, 
we assume that the electronic cloud in the H — O bonds have the same eigenfrequency uj p 
as the H — O ensemble. In the static limit x = q e D/m e ujp (polarizability a = ql/m e uip 
in p = aD), and we get a polarization P = p/a^ = qlD/m e aQLUp, where ao is of the 
order of the atomic size. We get an internal field E int — —A-kP — — {Airq^/m e aQ) D/lo^ = 
— (cjg/cjp) D, where u>o is a frequency of the order of atomic frequencies. Consequently, 
the dielectric function e in equation D = eE = e(D + Ei nt ) is given by e ~ — oj^/oSq 
(ujp/ujQ <C 1), which is precisely the static dielectric function given by equation (fl9| . 
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The dielectric properties of water are still a matter of debate. It is agreed 
that the permitivity dispersion of water is described to some extent by a 
Debye model of the form e = a + 6/(1 — iuor), where a and b are semi- 
empirical parameters and r ~ r/aVT is a relaxation time; r\ denotes the 
viscosity and T is the temperature^ This Debye model assumes mainly an 
orientational polarizability of electric dipoles, which, due to the preservation 
of the directional character of the O — H bonds, is compatible with the 
plasma model suggested here for water. Therefore, the contribution given by 
equation f fT9l ) should be added to the above Debye formula for the dielectric 
function, which becomes 



u> 2 -u> 2 



e = a + — + — — -| . (20) 

1 — IUJT UJ 2 + UJq 

Parameters a and b in equation ([20]) are related to the static permitivity e 
and high-frequency permitivity £ m through 

e = a + b - ujI/ujI , Boo = a + 1 . (21) 

We may neglect uj 2 /uq here because it is too small, and we may also take 
eoo = l(a = 0). The static permitivity Eq = b is given mainly by the electric 
dipoles. Let p be such an electric dipole. Its energy in an electric field D 
is — pD cos 9, where 9 is the angle between p and D. The thermal distribu- 
tion of such dipoles is dw ~ exp(— pD cos 9 /T)d(cos 9), .where T denotes the 
temperature. We get easily the thermal average (cos^) = —L(pD/T), where 
L(x) = cothx — 1/x is the well-known Langevin's function. 

We take p = 2ez e (a/2) = ez e a, where a ~ 3A and delocalized reduced 

charge associated with the H — O dipole. We estimate the argument pD/T of 

the Langevin's function. At room temperature, we find pD/T ~ 3 x 10~ 4 Dz e . 

For pD IT = 1 this corresponds to an external field D = ^- x 10 4 esw, or 
1 — 1 

D = 10 8 / z e V /mllll This is an extremely high field, so we are justified to take 
pD/T 1, and L{pD/T) ~ pD/3T. We get therefore a polarization P = 
—np (cos 9) = np 2 D/3T, an internal field E int = —A-kP = —4:iinp 2 D/3T, and 
a permitivity 

g ° = 6= l-W/3T (22) 

from D = eE = e(D + Ei nt ). This is the well-known Kirkwood formulal^l 
For the empirical value eo = 80, we get (at room temperature) a reduced 



10 See, for instance, H. Frohlich, Theory of Dielectrics, Oxford (1958); P. Debye, Polar 
Molecules, Dover, NY (1945). 

11 lesit = 3 x 10 4 V/m, J. D. Jackson, Classical Electrodynamics, Wiley, NJ (1999). 
12 See, for instance, H. Frohlich, loc cit. 



7 



charge z e ~ 10 2 . This is in good agreement with the H +z — O 2z plasma 
charge z estimated above. 

Cohesion and thermodynamics. Recently, a model of liquid has been 
introduced^ based on an excitation spectrum (per particle) of the form e n = 
—e + £i(n + 1/2), where e is a cohesion energy and e\ is the quanta of 
energy of a harmonic oscillator with one degree of freedom; n represents here 
the quantum number. The model includes also the kinematic correlations 
(spatial restrictions) of the movement of the liquid molecules. This model 
leads to a consistent thermodynamics for liquids, arising from a statistics 
which is equivalent with the statistics of bosons in two dimensions. 

For water, the cohesion energy per particle Eq can be estimated from the 
vaporization heat (~ 40A; J/mol). It gives eo ~ 10 3 K. On the other hand, it 
was shown in a previous papeo that the transition temperature beween a 
gas and a liquid of identical particles is approximately given by 

Tt = tnSln) ' (23) 

where T = h 2 n 2 ^ 3 /m is a gas characteristic temperature. We can apply this 
formula to water disssociation, taking n as the density of hydrogen atoms, 
m as the mass of two hydrogen atoms and T t = 383K (at normal pressure; 
eo depends on the inter-particle spacing). We may neglect the oxygen, as 
it is too heavy in comparison with the hydrogen atoms. We get T ~ 2K 
and the above formula gives e ~ 2000.fr ~ 200meV for the cohesion energy 
of water per molecule, which is consistent with the above estimate (leV ~ 
11.6 x 10 3 if; n ~ 1/a 3 with a = 3A and h ~ 10~ 27 erg ■ s; Bohr radius 
an = h 2 /m e e 2 ~ 0.53A, e 2 /an — 27.2eV, where m e is the electron mass)0 

The plasma oscillations obtained above can be quantized and the energy 
levels of the plasma read 

E n = J2 + 1/2) = 7^3 y?c ■ ^v( n + V2) , (24) 

q 

where q c is a cutoff wavevector. The prefactor in equation (i24l ) is Vq 3 /6n 2 ~ 
iV(ag c /4) 3 , so the energy levels given above can be written as 

E n = Ne l (n + 1/2) , (25) 



13 M. Apostol, J.Theor. Phys. 125 163 (2006). 

14 M. Apostol, Mod. Phys. Let. B21 893 (2007); see also M. Apostol, J. Theor. Phys. 
123 155 (2006). 

15 It is worth noting that the mechanism of vaporization assumed here implies the dis- 
sociation of the water molecule. 



8 



where e\ = (aq c /4)huj p . These energy levels correspond to a harmonic os- 
cillator with one degree of freeedom. It follows that the present description 
of water as a two-species of highly dissociated ionic plasma provides a fur- 
ther support for the liquid model mentioned above. If we take q c ~ 1/a the 
energy quanta E\ = (ag c /4) 3 /iu;p =~ 3zmeV represents the E\ parameter in 
the spectrum of the liquid. (The plasma frequency given by equation (TTT1) is 
uj p ~ 200zmeV). 

Debye screening and the correlation energy. As it is well-known the 
plasma excitations described above represent collective oscillations of the 
density in the long wavelength limit. At the same time they induce correla- 
tions in the ionic movements. For a classical plasma these correlations are 
associated with a screening length given by the Debye-Huckel theory as 

k- 1 = (T/2Annz 2 e 2 ) 1/2 (26) 

for our case (k^ 1 = (T /Arce 2 Ylii n i z f) 1 where i labels the ionic species with 
density rii and charge ezj). The formula is valid for the Coulomb energy 
z 2 e 2 /a much lower than the temperature T. In the present case we have 
z 2 e 2 /a ~ 45K (for z ~ 3 x 1CT 2 ), which shows that the above condition is 
fulfilled. From (1261) we get k~ x ~ lA (at room temperature), in agreement 
with the present molecular-dissociation model. The correlation energy per 
particle is given by 




{tcorr = — (e 2 / a) a/ 7i"e 2 /Ta(^ i riiZ 2 ) 3 ^ 2 ) . The estimation of this energy gives 
£ C orr ~ 10 2 K (at room temperature). It contributes to the cohesion energy. 

Sound anomaly. The sound-like branch oj 2 — oj s = v s q, where v s = 
^/9nx/ (M + 2m) according to equation (TT2l . is distinct from the ordinary 
hydrodynamic sound whose velocity is given by the well-known formula 
t> = 1/ y/nnm for a one-component fluid, where n is the adiabatic com- 
pressibility. For the present two-component fluid (H +z — 0~ 2z plasma), the 
velocity of the ordinary sound is given by v q = 1/ a/ ku(M + 2m). The former 
represents a non-equilibrium elementary excitation, whose velocity v s does 
not depend on temperature, while the latter proceeds by thermodynamic, 
equilibrium, adiabatic processes, and its velocity t>o depends on temperature 
thorugh the adiabatic conpresibility k. In order to distinguish them from the 
hydrodynamic sound we propose to call the sound-like excitations derived 

16 See, for instance, L. Landau and E. Lifshitz, Course of Theoretical Physics, vol. 5, 
Statistical Physics, Elsevier (1980). 
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here density "kinetic" modes or "densitons". The distinction between the 
two sounds is made by a threshold wavevector q t in the following manner. 
Suppose that there is a finite lifetime r for the sound-like excitations u> s prop- 
agating with a velocity v s and a corresponding meanfree path A = v s t. If 
the sound-like wavelength A is much longer than the meanfree path, A 3> A, 
then we are in the collision-like regime (uj s t <C 1), and the collisions may re- 
store the thermodynamic equilibrium. In this case the hydrodynamic sound 
propagates, and the sound-like excitations do not. This condition defines the 
threshold wavevector q t = l/v s r. In the opposite case, q^> q t (collision-less 
regime), it is the sound-like excitations that propagate, and not the hydro- 
dynamic sound. The finite lifetime r originates in the residual interactions 
between the collective modes and the underlying motion of the individual 
particles. It is easy to estimate this residual interaction^ It is given by 
y/eT, where e is the mean energy per particle corresponding to the motion 
of the individual particles. We get therefore r ~ %j \[€F and the threshold 
wavevector q t = \feT/hv s . It is difficult to have a reliable estimation of the 
mean energy e; for a resonable value e = WmeV we get q t — O.lA at room 
temperature for v = 3000m/s, which is in good agreement with experimental 
data. 

Indeed, the phenomenon of two-sound anomaly in water is well-documented 1^1 
Neutron, X-ray, Brillouin or ultraviolet light scattering on water revealed the 
existence of a hydrodynamic sound propagating with velocity Vq ~ 1500m/s 
for smaller wavevectors and an additional sound propagating with velocity 
~ 3000m/s for larger wavevectors. In addition, though both sound veloc- 
ities do exhibit an isotopic effect, their ratio does not. According to the 
above discussion, we assign this additional, faster sound to the sound-like 
excitations derived here. We can see that both vq and v s given above ex- 
hibit a weak isotopic effect, while their ratio v s /vq = 3ny / K% does not. 
From v s = \J$nx/ (M + 2m) = 3000m/s we get the short-range interaction 

X — 7eV ■ A 3 . Similar results are obtained for other forms of dissociation of 
the water molecule, like OH- - H + or OH- - H 3 + , so the H +z - 0~ 2z 
plasma model employed here can be viewed as an average, effective model 
for various plasma components that may exist in water. 

Another possible anomalous sound. It is worth calculating the spectrum 
given by equations of motion ( fTOl ) without neglecting higher-order contribu- 

17 M. Apostol, Electron Liquid, apoma, MB (2000). 

18 See, for instance, J. Teixeira, M. C. Bellissent-Funel, S. H. Chen and B. B Dorner, 
Phys. Rev. Lett. 54 2681 (1985); S. C. Santucci, D. Fioretto, L. Comez, A. Gessini and 
C. Maschiovecchio, Phys. Rev. Lett. 97 225701 (2006) and references therein. 
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Figure 1: The spectrum of the density oscillations given by equation ( 1251 ) for 
the H +z — 0~ 2z plasma with the same short-range interaction between ionic 
species. 



tions in q 2 . The result of this calculation is given by 



2 2 

^7 > = -u. 



1 + Ax 2 ± Vl + 2Bx 2 + A 2 x 4 



(28) 



where 



A = -!-(2 + 5a + 2a 2 ) , B = ^-(2 - 13a + 2a 2 ) , a = m/M (29) 
9a 9a 



and x = v s q/u) p . It is shown in Fig. 1. 

Frequency uj 2 in equation (1251) represents the sound-like branch, which goes 
like u>2 — oj s = v s q in the long wavelength limit and approaches the horizontal 
asymptote uj 2 = lu p / VA ~ uj p ^Jm/2M for shorter wavelengths. Frequency 
uj\ in equation (1251) represents the plasmonic branch {ui ~ u p for q — > 0). In 
the long wavelength limit it goes like 

(M — m) 2 

Ui ~ u p + 9mM ^ g2// ^ ' g ~" ° • ( 3 °) 

Due to the large disparity between the two masses m and M we can see that 
the plasma frequency has an abrupt increase toward the short-wavelength 
oblique asymptote given by 

u a ~ \fAv s q ~ y/2M/9m + 5/9v s q . (31) 
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For small values of u p (vanishing Coulomb coupling, z — > 0) this asymptotic 
frequency may look like an anomalous sound propagating with velocity 



v a ~ ^2M/9m + 5/9v s . (32) 

For water, we get v a ~ 2v s from this formula. However, the ratios v a /v s 
or v a /vo exhibit a rather strong isotopic effect, which is not supported by 
experimental data. 

Multi-component plasma. The model presented herein might be general- 
ized to a multi-component plasma consisting of several ionic species labelled 
by i, each with number iVj of particles, density n,, charge z^e and mass m^, 
such that Ylii z i n i = 0- 

The lagrangian of the density oscillations is given by 

L = ~h E iq m i^^(q)^(-q) + ^ V,,,, n,n,q 2 [vM + x(q)] «i(q)«j(-q)+ 

+^Eiq n i^#(q)«i(-q) > 

(33) 

where <Pij(q) = inZiZje 2 /q 2 . The equations of motion are given by 



rriiU 



Hi + 4:7ie 2 Zi z j n j u j + q 2 X n j u j — —iqezi(p . (34) 



Making use of the notations 

Si = S ^2z 2 n i /m i , S 2 = S ^n i /m i , S 3 = y^ZjUi/mi , (35) 



the eigenfrequencies of the system of equations ( 13411 in the long wave- 
length limit are given by 

^^ = WS 1 = £^^ , (36) 

which represents the plasma branch of the spectrum, and 

u 2 2 ~ co 2 = (S 2 - S 2 3 /Si) X q 2 = v 2 s q 2 , (37) 

which represents the sound-like excitations^ The plasma branch of the 
spectrum has an oblique asimptote given by U\ ~ uo a = y/x^2Qj which may 



19 The sound velocity given by |37|) is always a real quantity, as a consequence of the 
Schwarz-Cauchy inequality. 
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be taken as an anomalous sound propagating with velocity v a = \fxS~2 for 
small values of ui p . The ratio of the two sound velocities is given by 

1 , . 

Va/Vs = , (38) 



which is always higher than unity. The sound branch of the spectrum has 
an horizontal asymptote given by U2 — \A — Sf/SiS^Up. For the H +z — 
0~ 2z plasma we can check from ( 1381) that v a /v s ~ (2M/9m + 5/9) 1 / 2 ~ 2, 
and u 2 — 3^/m/2Mw p , as obtained above. As we have discussed above 
this ratio exhibits a rather strong isotopic effect, which is not in accord 
with experimental data. We assign therefore the additional sound to sound- 
like excitations propagating with velocity v s given by equation (1371) . The 
ordinary, hydrodynamic sound in a multi-component mixture has the velocity 



Vq = 1/ y/ k riirrii. It can be shown that v 2 /v 2 > n 2 nx for a neutral multi- 
component mixture. 

The internal field is given by 

E int = -47re ; (39) 

i 

we get easily from equations (1341 ) 



E int = -iq<P 2 Up 2 (40) 

and the dielectric function e = 1 — uo 2 /uj 2 , as expected. 
Structure factor. The structure factor is defined by 

S(q,u) = ±J dvdr'dt (5n{r, t)5n{r', 0)) e ^^')-^ = 

(41) 

- N fdt(6n(q,t)6n(-q,0))e- i(Jt , 



2im 



where the brackets stand for the thermal average (we leave aside the central 
peak). Since 

5n(q, t) = -iq y~] n^(q, t) , (42) 



it becomes 



S(q, u) = ^ [ dt J2 Uin i M*) u i(°)> e ~ iUt > ( 43 ) 

ij 
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where we dropped out the argument q. 

In order to calculate the thermal averages we turn back to the system of 
equations ( 1341 ) without the external electric field. This system can be written 

as 

(-lu 2 + aS^x + bS 3 y = , 

(44) 

aS 3 x + (-u 2 + bS 2 )y = , 
where a = 4-ne 2 , b = xg 2 , 51,2,3 are given by equation ([35]) and 



x 



n ^— ' 



y 



n ^— ' 



riiUi 



(45) 



In addition, 



anzi bn 
Ui = yX H -y . (46) 

The system of equations (14"4"|) has two eigenfrequencies Ui )2 as given by equa- 
tions (l36l) and (j3*7l) . The corresponding eigenvectors are given by 

(47) 



x\ ~ 5*1 , yi ~ 53 ; x 2 ~ 65 3 , y 2 ~ — a^i 



in the long wavelength limit. According to equation ( 1461 ) the coordinates m 
can be written as 



u 



(1,2) 



anzj 



rrijU 



"3^1,26 + 



6n 



1,2 



— —y^ e 

rriiuU 



(48) 



and one can see that they are coordinates of linear harmonic oscillators 

/2 . The ther- 



u 



(1.2) 



with frequencies 0*1,2 and potential energies mi(jj 2 2 

mal distribution of the coordinate u for such an oscillator is given by dw 
^Jmuj 2 /2nT exp {—muo 2 u 2 /2T) du in the classsical limit, where T denotes the 
temperature (T 3> faj). It follows 



(1.2) (1,2) 



T 



2 



(49) 



Writing 



(1) _iu>it 



(2) ia<2t 



tii = u- 'e~— + ure^ (50) 

and making use of equation (1491) the structure factor given by equation (T43l 
becomes 



5(g,^) = iVTg 2 K> 2 / 



n 2 m,- 



(51) 
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We can see from this equation that the relevant sound contributions are given 
by 

S (Q^) - ~T~ (^2 n 2 i/n 2 mi \ 6(u - v s , a q) . (52) 

V s,a \ i J 

Asymmetric short-range interaction. Up to now, the short-range in- 
teraction was assumed to be the same for all ionic species. In general, 
we may introduce a short-range interaction Xij depending on the nature of 
the ionic species. If this interaction is separable, the solution given above 
for a multi-component plasma holds with minor modifications. For a non- 
separable short-range interaction, appreciable changes may appear in the 
spectrum, which may exhibit multiple branches. Such a spectrum may serve 
to identify the nature (mass, charge) of various molecular aggregates in a 
multi-component plasma. It is worth noting that a range of frequencies 
10 10 s -1 — 10 12 s _1 is documented in living cells by microwave, Raman and 
optical spectroscopies and by cell-biology studies, upon which the theory of 
coherence domains in living matter is built @ 

We consider here again the H +z — 0~ 2z plasma with different short-range 
interaction Xhh — Xi,Xoo — X2 , Xoh — X3? it still exhibits two branches of 
frequencies, a plasmonic one (ooi) and a sound-like one (002), but the spectrum 
may have certain peculiarities (the dielectric constant is not affected by this 
modification). Equations of motion (fT5l ) become now 

mil + 2nq 2 ((p + Xi) u ~ nq 2 (2ip — X3) v = —izeq<j) 

(53) 

Mv + nq 2 (4ip + X2) v — 2nq 2 (2ip — X3) u — 2izeq<j) . 
We introduce the notations 

a = 2nq 2 Lp/m = 87ine 2 z 2 /m , 61,2,3 = ^Xi, 2,3/^ • (54) 



The dispersion relations can be computed straightforwardly. In the long 
wavelength limit (q — > 0) we get the plasmonic branch 

ul~(l + 2a)a+ — — q , (55) 

1 + la. 

20 See, for instance, H. Frohlich, Phys. Lett. A26 402 (1968); Int. J. Quant. Chem. 2 
641 (1968); S. J. Webb, M. E. Stoneham and H. Frohlich, Phys.Lett A63 407 (1977); S. 
Webb, Phys. Reps. 60 201 (1980); S. Rowlands et al, Phys. Lett. A82 436 (1981); S. C. 
Roy, Phys. Lett. A83 142 (1981); E. del Giudice et al, Nucl. Phys. B275 185 (1986). 
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v s % /to p 



v q/(0 



Figure 2: Excitation spectrum given by equation (l59l) for the H +z — O 2z 
plasma with short-range potentials xoo — Xhh — and Xoh — X 7^ 0. 



where (1 + 2a)a 
branch 



IQime 2 z 2 / \i is the plasma frequency, and the sound-like 
a(46i + b 2 + 46 3 ) 2 



-9 



2 2 



1 + 2a 

one can see that the sound velocity v s is always a real quantity. 



(56) 



The sound-like branch exhibits an asymptote in the short-wavelength limit 
given by 

' 2bi + ab 2 - x /(26i - ab 2 ) 2 + 8ab 2 <r 



uj 2 ~ 



(57) 



whose slope may have either sign or vanish. It is easy to see that this slope 
is positive for b\ < b\b 2 , negative for b 2 > b\b 2 (when the sound-like branch 
has a maximum value) and it vanishes for 6| = b\b 2 (when the sound-like 
branch has an horizontal asymptote). In the case of a negative slope the 
sound velocity may exhibit a negative velocity and the sound may suffer a 
strong absorption for moderate values of the wavevector, which may indicate 
an anomalous or unphysical situation. 



We return now to the plasmon branch given by equation (155!) . and write it 
as 

2x 2 — AaXx + a 2 9 

-q , 



wl + b 2 - 



i + 2a - • (68) 

where A 2 = b\jb\b 2 and x = \Jbi/b 2 . It is easy to see that for A 2 > 1 the 
plasmonic spectrum exhibits a dip around a certain value go of the wavevector 
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q for ^A — a/A 2 — l/2j a < \/bi/b 2 < ^A + a/A 2 — 1/2 j a; it approaches an 
asymptote with a positive slope for q — ► oo, which may define again an 



anomalous sound for small values of u p . 



We illustrate these anomalies for a particular case of short-range interaction 
Xi,2 = and X3 = X (&3 — n xl m )- The dispersion relations of the system of 
equations f l53l) become 



2 1 2 

^1,2 = 



Cl 4- 2«) 2 

i ± \/ 1 - 4^v/^ 2 + 1 \ a ] *tf/<4 



(59) 



The plasmonic branch has a minimum value for go — 2^/m/Mwp/ii s , where 
the sound-like branch has a maximum value (~ \j2m/Mu) p ). The spectrum 
is shown in Fig. 2. Using u p ~ lO 13 ^ 1 estimated above and the sound 
velocity v s ~ 3000m/s in water we get g^ 1 ~ 6A. We may expand u)\ 
in series of (q — q ) 2 around its mimimum value at q and get uj\ ~ u p + 
(M/4m + l)(t»^ 2 /^)(g - g ) 2 = uo p + (1 + 4m/M)v 2 8 {q - g )7^V This is 
similar with the rotons-like dispersion relation discussed in connection with 
the coherence domains in water I 2 *! Although this might be an interesting 
suggestion, it is inconsequential here, because ui p is too small in comparison 
with the temperatures at which water exists and, therefore, this "dip" feature 
has no effect for the water thermodynamics. 

Conclusion. We summarize the main features of the model suggested here 
for liquid water. First, we assume, as it is generally accepted, the four, direc- 
tional sp 3 -oxygen electronic orbitals. The electron derealization along two 
such orbitals together with a corresponding derealization of the hydrogen 
electronic charge lead to the water cohesion. It is represented by the cohesion 
energy e discussed here. Within such a picture, we can still visualize the 
oxygen and the hydrogen as neutral atoms, moving around almost freely (as 
a consequence of the uniformity of the environment; this gives a noteworthy 
support to the "hydrogen bonds" conce P t)B To this picture the present 
model adds another component, arising from a very small charge transfer 



21 G. Preparata, QED Coherence in Matter, World Sci (1995). 

22 The point of view taken in this paper is that the hydrogen bonds in water are intro- 
duced in order to account for the uniformity of the environment of a water molecule in 
liquid water. As such, it helps understand the cohesion. However, a consistent upholding 
of the hydrgen-bonds concept would mean a vanishing dipole momentum of liqud water. 
Pauling himself, (L. Pauling, loc cit) who introduced originary this concept, qualifies it by 
admiting an asymmetry in the four hydrogen bonds around an oxygen ion, arising from 
the two-out-of-four occupied orbitals. We suggest that the uniformity of the environment 
makes the hydrogen atoms (ions) moving as independent entities, while the asymmetry 
induces a small charge z, so the ion motion is subjected to Coulomb (and short-range in- 
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between hydrogen and oxygen atoms, leading to a H +z — 0~ 2z plasma, with 
the reduced charge z. It may originate in the weak asymmetry of the two 
occupied sp 3 -oxygen electronic orbitals with respect to the other two unoc- 
cupied orbitals. Under these circumstances, the hydrogen and oxygen ions 
interact, both by Coulomb and short-range potentials. This interaction gives 
the plasma frequency and the sound-like excitations frequency. The plasmons 
contribute to the excitations which give rise to a consistent thermodynamics 
for liquids, in a model introduced recently. In addition, the ionic plasma os- 
cillations entail oscillations of the delocalized electronic cloud, with the same 
eigenfrequency. Subjected to an external field, these electronic oscillations 
produce an intrinsic polarizability which removes the u — singularity in the 
plasma dielectric function (the Uq frequency). In addition, the magnitude of 
the electric moment p which is responsible for the orientational, static di- 
electric function is in satisfactory agrement with the plasma charge z derived 
herein. 

On the basis of this model we are able to understand to some extent, both 
qualitatively and in some places even quantitatively, the sound anomaly, the 
dielectric function (permitivity dispersion), the structure factor, cohesion and 
thermodynamics of water. The model is extended to a multi-component clas- 
sical plasma, including an asymmetric short-range interaction between the 
components, which might be relevant for more complex structural aggregates 
like those in biological matter. 



teractions) . The electric moment is ascribed to the directional character of the sp 3 -oxygen 
electronic orbitals and the charge transfer between oxygen and hydrogen. Thereby, the 
hydrogen-bond concept is employed here through its two features, directionality and uni- 
formity, with a slight asymmetry, all viewed as independent qualitative ingredients. 
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